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^S) ' In this paper, we study the influence of quantum effects to chaotic dynam- 

o ■ 

\ ics, especially the influence of Pauli effect and dynamical symmetry breaking 

\ to chaotic motions. We apply the semiquantal theory to the Sp(6) fermion 



symmetry model in nuclear collective motion. We demonstrate that quantum 
chaos appears when dynamical symmetry is broken. We further show that 
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\ dynamical Pauli effect enhances quantum chaos. 
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I. INTRODUCTION 

Quantum chaos is one of the most challenge research topics in both theoretical and 
experimental physics. In the past decade, in order to understand the concept of quantum 
chaos, people have been focusing on seeking the generic behavior in quantum spectra and 
wave functions of various dynamical systems whose classical dynamics are chaotic. In these 
investigations, two significant features have been found. One of them is the relation between 
the quantum energy level statistical distributions of the classical chaotic systems and the 
GOE/GUE spectral statistical universality classes of the random matrix theory The 



other feature is the existence of "scars" in individual eigenstate wave functions, i.e., the 
appearance of strong locahzation of wave functions along the unstable periodic orbits in 
chaotic systems 0. However, the dynamical mechanism of quantum manifestation of chaos 
is still not well understood so far. 

In order to study the dynamical mechanism of quantum manifestation of chaos, one 
of us developed an approach, called the semiquantal approach P|. The main motivation 
of this approach is to provide a framework to study how quantum fluctuations, quantum 
correlations and other quantum phenomena manifest themselves in the classical trajectories. 
This approach differs from the above mentioned investigations. It starts from quantum 
systems without reference to classical limit and it directly explores the dynamics of quantum 
effects in classical chaotic trajectories. Therefore, it provides a systematic way to explore 
the intrinsic mechanism of quantum manifestation of chaos 

The general theoretical framework of the semiquantal approach is the stationary phase 
approximation of the generalized coherent state path integral |0. Each set of generalized 
coherent states associated with the dynamical group G of a quantum system forms an over- 
complete set of states for an irreducible representation space of G. The coherent states are 
also in one-to-one corresponding to a geometrical space which has a symplectic structure. We 
called this geometrical space the quantum phase space. Based on path integral formulation, 
coherent states then embed each quantum system in a quantum phase space, regardless of 
whether the system has a classical counterpart or not. The stationary phase approximation 
further results in a set of dynamical equations of motion on this quantum phase space, which 
are similar to the classical equations of motion. But the Hamiltonian function governed these 
dynamical equations of motion is the expectation values of the Hamiltonian operator in the 
generalized coherent states, which contain the leading quantum fluctuations and correlations. 
We called these equations of motion the semiquantal dynamical equations. This formulation 
allows us to directly explore the dynamical effects of quantum fluctuations and correlations 
to chaos 10]. 

Based on this formulation of semiquantal dynamics, we will study in this paper the 
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dynamical effect of the Pauli exclusive principle to chaotic motions. We shall take the Sp(6) 
model, the Sp(6) symmetry of the fermion dynamical symmetry model in nuclear many-body 
systems |]1T|, as our example. This model is a microscopic description of nuclear collective 
dynamics with fermions as the building blocks. A significant feature of this model is that the 
Pauli principle of fermions has a crucial influence on the low-lying nuclear collective modes, 
which manifests themselves in experimental data [jl2l. When we define a set of canonical 



deformation variables (the quantum phase space) associated with the Sp(6) coherent states. 



the connection between Bohr-Mottelson geometrical model and the Sp(6) model can be 



built [|14|. The Pauli effect in the Sp(6) model itself manifests in the nuclear deformation 



in terms of some constraints on the nuclear surface (deformation) variables [T^. Within 
the semiquantal framework, the Poincare's picture of chaos can be applied to describe the 
possible chaotic dynamics in nuclear collective motion in the aforementioned quantum phase 
space. Therefore, dynamical Pauli effect to chaos can be explicitly studied. Our result 
shows that dynamical Pauli effect enhances the pairing correlations (quantum correlation) 
and inhibits nuclear deformation. This gives rise to the breaking of dynamical symmetry 
and therefore leads to the occurrence of chaos. 

The paper is organized as follows: In Sec. 2, the nuclear shape in terms of the microscopic 
fermion degrees of freedom is parameterized through the Sp(6) coherent states. The coherent 
state of the Sp(6) model in the nuclear body frame is constructed so that the quantum 
phase space of the Sp(6) model is connected with the nuclear geometrical model. In Sec. 3, 
the semiquantal dynamical equations of motion for nuclear deformation is derived. With 
the consideration of time-reversal invariance and valence nucleon number conservation, the 
resulting equations of motion are reduced to a set of Hamilton-like equations in a four- 
dimensional quantum phase space, which describes the nuclear surface motions. In Sec. 4, 
we explore the geometrical manifestation of dynamical Pauli effect and numerically examine 
the dynamical effect of the Pauli principle and dynamical symmetry breaking in semiquantal 
framework, from which we find that quantum chaos appears when dynamical symmetry is 
broken and dynamical Pauli effect enhances quantum chaos. Finally, a summary is presented 



in Sec. 5. 



II. PHASE SPACE PARAMETERIZATION OF NUCLEAR SHAPE AND 
GEOMETRICAL INTERPRETATION 



We begin with the fermion dynamical symmetry Sp(6) model |TT| which describe the 
nuclear collective motions. The Sp(6) group is spanned by the collective pairing opera- 
tors {S, D^} and the collective particle-hole quadrupole operators {P\^] A = 0, 1, 2; = 
—A, A}. The model Hamiltonian operator is given by 

H = eh + GoS+S + biPi ■ Pi , (1) 

1=1,2 

where n is the particle number operator. The coefficient e, Gq and bi represent the strength 
of the internal energy, effective pairing and multiple interactions, respectively. It is worth 
noticing that the D-pairing interactions are not included in Eq. (|1|) because such interactions 
can be replaced by the S-pairing and multiple interactions in the Sp(6) symmetry Also, 



in order to study chaotic dynamics of nuclear deformation under the Sp(6) model, we should 
only consider the dynamics of different intrinsic bands rather then the collective motion 
within a single band. Therefore, we can let bi = 62 in Eq. (|l]). 

To study the nuclear surface motions in the Sp(6) model within the semiquantal frame- 
work, one must ffist construct the corresponding quantum phase space. It has been shown 
that such a quantum phase space is uniquely realized by the coset space USp(6)/U(3) 
through the Sp(6) coherent states. (The detailed construction of the coherent states of the 
Sp(6) model can be found in Ref. [jl4[. It has also been roughly summarized in Appendix 



A.) 

In the following subsections, we will parametrize the nuclear shape in the above men- 
tioned quantum phase space by visualizing nuclear deformation to be a geometrical descrip- 
tion of the microscopic Sp(6) fermion model. As it is well-known, the dynamics of nuclear 
deformation has time reversal invariance. Meantime, it is also convenient to describe nu- 
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clear deformation in the nuclear intrinsic frame in which adequate deformed variables can 
be chosen as the coordinates of the Sp(6) quantum phase space. 



A. The time reversal invariant deformed coherent states 

First, we shall look at the connection of the Sp(6) Coherent States (CS) with nuclear 
deformation. In the macroscopic picture of the nuclear geometrical model, the rotational 
and vibrational energy spectrum of nuclei are described in according to the nuclear surface 
motions. The nuclear surface can be defined by 



R{t, e,(f)) = Ro( aUt)Yx^,iG, (jy)) + Rc 

\ A=0,l,2,/i=-A,- -,A / 



(2) 



where Rq is the maximum thickness of a nucleus, Rc the radius of frozen core, and Yx^{0, 0) 
the spherical harmonic function. In Eq. (|^), the quadrupole terms (with A = 2) play the 
dominant role in the description of nuclear deformation. They determine the strength of 
nuclear surface motions. The monopole term (A = 0) describes the shell thickness of the 
undeformed compressible nucleus. When the center of nuclei is fixed in a reference frame 
(no translation), the effect of A = 1 term can be ignored. Besides, all the surface variables 
in Eq. (|^) satisfy the time reversal symmetry so that 

al^it) = {-ira,_,it) . (3) 

To generalize this concept into the framework of a microscopic description, we introduce 
a deformation matrix p, 

P = - |^aoo(t)l3x3+2 E^<(t)P2Mj , (4) 

which is related to the quadrupole momentum in the Sp(6) model, where Isxs is a 3x3 
identity matrix; and the submatrix block of the quadrupole momentum operator in 
the 6x6 matrix representation of the Sp(6) group. The detailed expressions can be found 
in Appendix A of Ref. [jl4|. These submatrices P^^ form an irreducible representation of the 
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unitary subgroup U(3). They determine the deformation properties of nucleus. The surface 
variables still obey the time reversal invariance condition (|), and the coefficient aoo(^) 
is chosen to be negative. This is because the shell thickness of the undeformed nucleus 
is positive, and also the quadrupole momentum must also be positive when the particle 
number is less than f2, where Vt is the pair degeneracy of normal-parity levels in a major 
shell (i.e. VLi in the Sp(6) model). In this paper we only consider the case of valence particle 
number n < Vt. For the case n > fi, an analogous discussion can be made because of the 
particle-hole symmetry. 

Next, we shall build the connection between the surface variables a\^{t) and the param- 
eters in the Sp(6) coherent states. As one can find in Ref. [1^ (also see Appendix A), the 
Sp(6) coherent state |?7) is completely determined by a 3 x 3 parameter matrix X. Meantime, 
the 3x3 hermitian matrix vXX+ can be spanned by 'P2^l'■ 

VXX^= E (3Ux,>^,it))P,,, (5) 

A=0,2,/t=-A,..,A ^ 

where Px^ obey the condition: Px-^ = (— 1)^P^^. Furthermore, the hermitian property of 
the matrix VXX+ implies that 

(^1-, = i-lTPx, (6) 

which corresponds to the time reversal invariance. 

Without any loss of generality, we can set the matrix -\/XX+ to be equal to the de- 
formation matrix p by letting Px^if) = {—l)'^ax^{t). The connection between the surface 
variables axij.{t) and the coherent state parameters in \t]) can then be carried out. The 
properties of nuclear deformation determined by the Sp(6) CS are described by Eqs. 
and the requirement vXX^ = p. The explicit form of the matrix vXX^ is given as 



Vxx+ = p = - 



"00 + 



|«20 



2 "21 



f"21 



"00 



- V2> 



"20 



-jf"21 



^3, 



a 



22 



2 "21 



(7) 



"do + 



|"20 



6 



B. The deformed coherent states in the intrinsic frame (ICS) 



In order to make the dynamics of nuclear deformation manifestation, it is convenient to 
express the nuclear deformation parameters in the nuclear intrinsic frame (i.e. the nuclear 
body frame). We call the CS describing the nuclear deformation in the nuclear intrinsic 
frame the intrinsic coherent states (ICS). In the intrinsic frame, the deformation matrix is 
positive and diagonal. To deduce the ICS, a basic idea is to decompose the previous spherical 
tensor operators into proper cartesian tensor operators. 

The cartesian tensor operator dij are defined as follows: 



E {UiiBimU^j) , i,j = 1,2,3, 



(8) 



lm=l,2;S 

where the index 1 x,2 ^ z,3 y for dij. Such an index assignment corresponds to a 
choice of the left-hand coordinate system in order to be consistent with the notations used 



in Ref. [T^. The operators Bi^ are uncoupled particle-hole operators in the Sp(6) model. 



The connection between Bi^ and can be found in Ref. |TJ] . The transformation matrix 
U is defined by: 



V2 



1 -1 

V2 . (9) 

-I -i 

With the above cartesian tensor operators, the rank-one spherical tensor operators can be 
expressed in terms of dij , 

\/3 

Pu ^ (ydzx dxz ~l~ i dzy i dyz) , 



Pw = i\ 77 (a 



yx "'xy j } 



(10) 



and the rank-two spherical tensor operators -P2^ (quadrupole momentum operators) , 
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„ ( ^xx 2dzz ~\~ d, 



xx ^'-^zz I '-^yy I 1 



P21 

P22 
P2-2 



^ V^zi; ~r ^xz ^ ^yz ^ ^zy) i 
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~^xx ~l~ ~l~ ^ ^3^2/ ~^ ^ ^yx) • 



(11) 



Let us now consider the matrix representation in terms of dij as the basis, 



dij — > Rep(dij) — e 



E,, 



-E. 



(12) 



where E.ij is a 3 x 3 matrix with +1 in the i-th row and the j-th column, and zero everywhere 
else. The deformation matrix p in this cartesian basis is then expressed as 



Q'oo + \l\cx.20 

^J\{o^2l + q;2i) aoo - V^0i2o -i^{oi2i - 0:21) 



(13) 



f (^22 - a22) Vtv«21 - «2l) 

+ V!(«22 + a^2) 

where the parameters {ctA/i} maintain the tensorial structure under rotational transforma- 
tion, i.e. 



R \uj) p {ax^,) R{uj) = p (a'^^) , 



(14) 



and a'x^ is given by 



(U=— A 



(15) 
(16) 



with R{u;) being a Euler rotational matrix, u the Euler angular, and D^^{ijj) the Wigner-£) 
function. One can always find an Euler rotation R(uj) such that the nucleus is described in 
an intrinsic frame in which p is diagonal (i.e. a'21 — and q;22 real): 



\S^\ 
\S,\ 

\s„ 



(17) 



Since pm = y^inX.in, the matrix Xj„ (the Sp(6) CS parameter matrix in the intrinsic frame) 
can be chosen as: 





S, 
-S„ 



(18) 



The minus sign of Sy in the above equation is due to the use of the right hand coordinate 
system. The exphcit form of the parameters Si is given as follows, 



1 



- VSa22j e 



Sx = - I aoo + "^"20 



ifx {aoO ,"20 ,"22 ;7i'00,'r20 ,''"22) 



So, 



«00 



V2 



^20 



(19) 



where the detailed functions of f^, fz and L can be found in Appendix C. In the intrinsic 



frame, the deformation degrees of freedom (vXX^) are reduced to three because the intrin- 
sic frame is chosen to be along the three principal axes of the nucleus. Together with the 
corresponding conjugate momentum variables, the phase space of nuclear deformation is a 
six-dimensional space. The six independent variables are now completely determined by the 
three complex parameters of X^ in ICS. 

Based on ICS, the nuclear deformation is described by the variables l^il only. The 
physical picture of the variables \Si\ can be seen from Eq. (0) which can be reexpressed in 
terms of cartesian coordinates as follows: 



R{t, X, z, y) = Rc + -Ro(|aoo| + 



xz xy 

4- Ofy I- 2n' -— -\-2n' 

I ^^xz ^2 ^ ^^xy ^2 ' ^^yz ^2 J ' 



(20) 
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where r = \/x'^ + z'^ + y'^. Since the intrinsic frame is chosen to be the three principal axes 
of nucleus in the intrinsic frame, the terms ct^^j;, and ayz vanish. The deformation matrix 
then becomes: 

Pin = -aoo(t)hx3 + (^ii(t)-^ii, (21) 

i=x,y,z 

where An is the subblock matrix in the matrix representation of dij, and 

anit) = aooit) + \Si{t)\ . (22) 

According to Eq. (pi]), the length of the nucleus along the ith-axis is given by Ro\Si(t) \ +Rc. 
In other words, {Si} are the variables measuring nuclear deformation. These are exactly the 
dynamical variables we want to use to study the nuclear surface motions. 

Once we have constructed the Sp(6) CS in the intrinsic frame, we can now derive the 
semiquantal dynamical equations of the nuclear collective motions. 



III. THE SEMIQUANTAL DYNAMICS FOR SP(6) MODEL 

The semiquantal dynamics of the Sp(6) model is determined by the time-dependent 
variational approach to the effective action (i.e. the stationary phase approximation of the 
generalized CS path integral): 

S = J dt{lCS\idt - H\ICS) . (23) 

The result of such a variation is a set of Hamilton-like equations 

^n^s* - 

'^3^' ~ as, ' 

^^^-^i = ' i = x,y,z. (24) 

These equations of motion are the semiquantal dynamical equations which approximately 
describe the quantum dynamics of nuclear deformation in the Sp(6) model. The word 
"approximately" (or semiquantal) means that this is an approximation of the Schodinger 
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equation described by the Hamilton-like equations of motion but it is derived purely from 
quantum mechanics. And the Hamiltonian function {H) = (ICS | if | ICS) contains all the 
leading order quantum correlations. In nuclear many-body systems, this set of equations of 
motion is indeed the familiar time-dependent Hartree-Fock-Bogoliubov equations under the 
constraint of dynamical symmetry group Sp(6) [ P^ . 

Note that Eq. (p4|) contains six equations of motion. However, not all these six equations 
of motions are independent. This is because for a given nucleus, the valence nucleon number 
is fixed: 

n={h) = ^ {\S,\^ + \Sy\^ + = 2n {alo + al^ + 2al^) . (25) 

In order to solve Eq. (p4D, the constraint of the fixing valence nucleon number must be 
imposed. In the next subsection, we shall first calculate the Hamiltonian function, and then 
introduce several canonical transformations to eliminate this constraint. 

A. The Hamiltonian function of the Sp(6) model 

The expectation value of the Hamiltonian operator ([l|) in the ICS is defined by 

{H) = e{n) + GoiS+S) + 62 E i^i ' Pi) ■ (26) 

1=1,2 

The detailed form will be given later. Since the valence nucleon number is conservative, the 
first term in Eq. (^) can be removed. Hence, Eq. (pB]) can be simplified as 

(H) = £'[su(2)] + -E'[su(3)] , (27) 

with 

-E'[su(2)] = Go{S'^S) , 

^[su(3)] = &2 E (Pi ■ Pi) ' (28) 

1=1,2 

where the coherent state expectation values can be exphcitly calculated in ICS, 
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(A ■ A) 
(A ■ A) 



and 



3n f2 

Kr, Al , 

7n Q^^ Q^^ (2 

+ —Ki + —7^2 - 

4 4^ 2 ^ 3 V 3 



\ n (I 



^2 = EI'^> 



1 ISiPA/l liSjpiSjSj 



4 

1 1 ; 



(29) 



(30) 



In Appendix B, we have shown that the Sp(6) model in this parameter representation (quan- 
tum phase space) contains all the general properties of the nuclear geometrical model [1^ . 



B. The semiquantal dynamical equations without constraint in the Sp(6) model 



To remove the constraint of the fixed valence nucleon number, several canonical trans- 
formations have to be made. Let the complex variables Si = qi + ipi, where and Pi are 
real variables. The Hamilton-like equation (p^) can then be reexpressed as: 



2n 



d{H) 
dqi 



2VL . _ d{H) 
3 dpi 



t = x,y,z, 



(31) 



where the Planck constant h has been absorbed into the Hamiltonian function. Next, we 
introduce a canonical transformation to generate the new canonical variables {Vi, Qi}, and 
choose the second type generating functions as 

.(2) 



Fr= E (/.(+) or /.(-)) 

i=x,y,z 



/.( + ) 



Qi 



sm 



V, 2 

Si 
V, 



(32) 
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This type of generating functions has eight assembles corresponding to positive or negative 

(2) 

Pi, i = x,y,z. To be exphcit, let us consider the special case, F| = J2i /«(+)• From this 
canonical transformation, we can find the connections between the new and odd canonical 
variables. 



a-^-^.f--^-?). (33) 



or 



dq, 



q- = -p. cos — 

' i 



p,^\V,sm^\. (34) 
This transformation is valid only for positive pi. In another case 

= J2i can get 



the same transformation but for negative pi. This is because in this case Pi = —yVf — qf = 
— \Vi sin -p^ I and = Vi cos -p^. All the eight generating functions transform the old variables 
{Px, Py, Pz) = {+,+,+) or (+,+,-) or (-,-,+) ... into new variables (V^, Vy, V^), where 
(+,+,+) means that px, Py and pz are all positive. Their connections to the variables Si 
are given by 

^+P^ = \S.r = rf, (35) 

and 

Si = qi + ipi = Vi exp f-^-^) . (36) 



V, 

Note that Eq. (|36D implies Vi ^ (if Qi ^ 0). In other words, the trajectories (Vx,Vy,Vz) 
cannot go through the section Vi = 0. Therefore we have eight types of independent initial 
conditions, namely, {Vx,Vy,Vz) G {(±, ±,±)}. Next, we will make the second canonical 
transformation. It transforms {Vi,Qi) into {Pi,Qi) such that Qz becomes an ignorable 
coordinate. We choose the second type generating function as 
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i^f = + QyPy ± Q.^P^ - - . 



(37) 



The connections between new and old variables are 



^P 



p. -p. -p 



1 = x,y, 



Vl = - Pi - Pi or P? 



El 12 _ Sn 



(38) 



By using these two steps of transformations, the relations between two sets of canonical 
variables {Si, S*} and {Pj, Qi} are given by 



\Sz\' 



p2 p2 

2Vt ^ ^ 



SxSy — PxPy GXp 



Qx Qy 
^Px Py / } 



The Hamiltonian function H can be expressed now in terms of {Px, Py,Qx, Qy)' 



where 



3n^ nVt 



9 



9 



(39) 



(40) 



with 



3n 



H.^^{pi+pi)-^[p:^+p;+pip^) , 

+ ^(l-P,^) (1-nW. COS (1^) 

+ V(l-^i) (l-nO^^^. COS - I) , 



20 y 



(41) 



(42) 



(43) 
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The new Hamilton-like equations become 

2n . _ dH 

— Qi = i^ — , i = x,y. (44) 

Note that the 3rd canonical momentum is a conservative quantity because Qz is an 
ignorable coordinate (i.e., it does not appear in the Hamiltonian function). Therefore, the 
quantum phase space of nuclear deformation is reduced to a four-dimensional space after 
solving the constraint of the fixing valence nucleon number. Thus, we have formulated the 
semiquantal dynamics of the Sp(6) model in a four- dimensional quantum phase space, which 
determines the dynamical evolution of nuclear deformation. 



IV. THE DYNAMICAL PAULI EFFECT AND CHAOTIC MOTION 



A. Pauli forbidden region in (Pj, Qj)— phase space 



Although we have formulated the semiquantal dynamics of the Sp(6) model in a four- 
dimensional phase space, this phase space is still an nontrivial geometrical space. As we have 
seen in the previous section, this quantum phase space is a compact space. It is embedded 
in a four- dimensional flat space. Besides, due to the Pauli effect of the Sp(6) model [|T^], not 
the whole domain of this quantum phase space is dynamically allowed. When the valence 



nucleon number n > some region is forbidden due to the Pauli effect |15|. Explicitly, we 



can see from Appendix A that the diagonal submatrix block parameters Y in the Sp(6) CS 
satisfy the condition Y^ = 1 — XX^. Consider the parameter space in the intrinsic frame. 
The matrix Y,-, 



±a/1 — Xj„Xj'„) has the form: 



/ 



± 



l-l^.l 









l-\Sz 









1^, 



(45) 



y\ ) 
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Since the matrix Yj„ is hermitian, it means that \Si\ < 1. Together with the condition of 
the fixing valence nucleon number, I]j = we obtain the following conclusion: When 
Q > n > ^ {oT J2k I'S'feP > 1) [see Eq. (p5|)], the condition l^il < 1 is not satisfied for the 
whole domain. A forbidden region appears in the quantum phase space, and its contour 
plotted in a2o — v^«22 plane can be determined by using Eq. (|B^) |T5[. When n < ^ 



(o^^ J2k I'S'fcP < 1), the condition l^il < 1 is always satisfied so that there is no forbidden 
region in the phase space. Note that the critical value n = ^ has an important consequence 
in the Sp(6) model. When n < the nuclear collective motions can be described in the 
symmetric representation of the subgroup SU(3) of Sp(6), which behaves the same as the 
bosonic SU(3) representation. However, when n > n > the Pauli principle forces the 
system to move to an asymmetric representation of SU(3). Correspondingly, the nuclear 
deformation will be limited in the quantum phase space. The condition \Si\ < 1 is indeed 
such a phase space constraint of the Pauli effect in the Sp(6) model 

In conclusion, within the constraint of the fixing valence nucleon number and the presence 
of the Pauli effect, the quantum phase space of the Sp(6) model is only a limited domain in 
a four- dimensional flat space restricted by |^ — P| — > and 1 — l^jp > 0. This is a pure 
quantum mechanical effect on geometry. In the following, we shall explore the dynamical 
consequence of these quantum effects to the chaotic motions. 



B. Analysis of chaotic motion in the Sp(6) model 

In this subsection, we shall numerically analyze the chaotic motion in the Sp(6) model. 
Since the semiquantal dynamics is determined by the Hamilton-like equations of motion 
(^), the chaos can be described in terms of Poincare' surface of section (PSS). We will flrst 
examine the nuclear surface motions in the dynamical symmetry limits of the Sp(6) model. 
The Hamiltonian in the intrinsic frame can be expressed as: 

1=1,2 

= Teffiaxf,, VTa^) + Veffiaxf,) 
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= T,ffiP.,Q,) + Vejf{P,), (46) 

where T^jf is an effective Icinetic energy which contains the general Icinetic energy and a 
velocity-dependent potential; and Vejf is an effective potential of the Sp(6) model derived 
in Appendix B. In the intrinsic frame, the Sp(6) model has two group chains: the SU(2) 
chain (with 62 = 0) and the SU(3) chain (with Gq = 0). These two limiting cases describe 
the nuclear vibrational mode and the 7-stable rotational mode, respectively. The numerical 
results of Fig. la and lb shows that their PSS are regular. This is because the existence 
of dynamical symmetries in these two limits guarantees the integrability of dynamics HlB 



However, when both Go and 62 are not equal to zero, the above dynamical symmetries are 
broken and chaos appears. This is shown in Fig. Ic. 

To demonstrate the chaotic behavior in the nuclear surface motions, we fix the rate Go/62 
and vary the particle number n to analyze the phase space structure of the semiquantal 
dynamics. There is a phase transition from the vibrational to the rotational modes in 
the nuclear surface motions. According to the result in Ref. [0, this phase transition 
is associated with the dynamical symmetry breaking from the SU(2) to the SU(3) limit 
with fixed G0/&2, and taking n as the control parameter. Therefore, chaos must occur 
accompanied with this phase transition [l^ll^. In the following numerical calculation, we 
take Gq = —.035 MeV and 62 = —-065 MeV with Q = 36. Note that the chaotic behavior 
also varies with the energy which can be regarded as a scale for analyzing the strength of 
chaotic motion. 

The numerical results show that the patterns of PSS in the Sp(6) model can be divided 
into three types with n increasing. Foe the cases of n = 2 ~ 8, the PSS have the 'similar' 
patterns with the energy increasing. We classify this type of PSS class A. The cases n = 
10 ~ 16 form class B, and = 18 ~ 36 form class C. Here, the word 'similar' means that the 
PSS are topologically similar, namely, the PSS of two states have the similar type of KAM 
surfaces and the similar distribution of chaotic phase space structure. For example, in Fig. 
2c and 2d, the state oi E = —6.8 MeV with n = 10 is similar to the state of E = —7.7 MeV 
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with n = 12. Analogously, the state E = —3.95 MeV with n = 10 is similar to E = —4 
MeV with n = 12, and E = -3.4 MeV with = 10 is similar to E = -3.5 MeV with 
n = 12. Physically, this classification depends on the contribution rate of the quadrupole 
and the pairing interactions. Roughly speaking, class A corresponds to the states with 
pairing interaction dominating, class C with quadrupole interaction dominating, and class 
B is the mixed states. In Fig. 2, six special cases of the PSS are plotted. They show the 
differences of PSS among the classes A, B and C. These Poincare surfaces of section (with 
Qx = 0) are plotted in the 020 — «22 plane (^a'22 = v^«22) , and have a 7 = | -axis symmetry. 

In order to define chaotic strength in a uniform way, we rescale the time t — ^ tVmax 
[Vmax = (maximum of |Ve//|)] such that the minimum values of the allowed energy are 
rescaled to —1 for all the values of n. This rescaling does not change the structure of the 
corresponding PSS. Also, for convenience, we introduce three energy scales: Eu, E^ and E^. 
The energy scale Es is the starting energy of the strong quantum chaos, Ey^ the starting 
energy of the weak quantum chaos, and E^ the upper bound of energy within which the 
Sp(6) model dominates the low- lying nuclear dynamics. Thus, for a given nucleus, the 
dynamical behavior of physical states can be divided into three types of motions: The pure 
regular motion which exists in the energy range from —1 to E^, the weak chaotic motion in 
the range of Ey, to Eg, and the strong chaotic motion in the range of Eg to Ey. Recalling 
the 'pattern classification' of the chaotic behavior aforementioned, systems with different n 
in the same class have similar patterns of PSS. Thus, Ey^, Eg and Ey are all n dependent. 
We plot the phase diagram of chaotic behavior in the nuclear surface motions in the E-n 
plane in Fig. 3. The chaotic strength is obtained by comparing the values of Eg and Ey^ in 
different n. 

From Fig. 3, we see that when n=2 and 4 (small n), the contribution of the quadrupole 
interaction (&2 Z]«=i,2(-P' ' -P') in Eq. (pB])) is much smaller than the pairing interaction 
{Go{S^S)). This implies that the system is near in the SU(2) C Sp(6) dynamical sym- 
metry. Correspondingly, the nuclear surface motion is near-integrable. Hence, for small 
n, the nuclear deformation is regular in all the energy ranges, as shown in Fig. 3 where 
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Ey^ c:^ Es — Eu for small n. This corresponds to class A. When n is near — , we find that 
the quadrupole interactions is dominant. The nuclei are approximately described by the 
SU(3) D Sp(6) dynamical symmetry. As we can see from Fig. 3, E^ and Eg are far from 
the minimum energy —1. Therefore the phase space is filled mostly by the regular motion 
as well (with n = 18 ~ 24). This corresponds to class C. When the particle number n is 
between 8 to 14, the pairing interaction and the quadrupole interactions compete against 
each other, and all the dynamical symmetries in the Sp(6) model are broken. This is the 
transition (mixed) region from the vibration to rotation in nuclear collective motions. The 
chaotic motion is very strong as it is shown in Fig. 3. This corresponds to class B. 

The most remarkable result we find in this study is in the region ^ < n < n. Naively, 
when n > the quadrupole interaction should be dominant so that the system should 
approximately have the SU(3) dynamical symmetry. However, the numerically result shown 
in Fig. 3 tells us that when n varies from ^ to the middle shell (^), the strength of chaotic 
dynamics is increased. This unusual result is indeed caused by the Pauli effect. The Pauli 
principle forces the system to jump from the symmetric SU(3) representation to the antisym- 
metric SU(3) representation when n > This reduces the nuclear deformation (inhibiting 
the quadrupole interaction and increasing pairing interaction). As a result, the system de- 
viates away from the SU(3) dynamical symmetry and therefore leads to the stronger chaos. 
This is the first example to dynamically examine the Pauli effect on chaos. It is also the 
first evidence of showing the enhancement of chaos in quantum system due to the Pauli 
principle. 



V. SUMMARY 

In this paper, the Sp(6) model have been chosen as an example to examine the chaotic 
dynamics in quantum systems caused by dynamical symmetry breaking and dynamical Pauli 
effect. It has been tested for many systems where dynamical symmetry breaking is closely 
related to quantum chaos P, p!6| , p!7| . However, this is the first paper to address the problem 
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of Pauli effect on chaotic dynamics. Since Pauli effect is a pure quantum medianical phe- 
nomenon and it has no classical limit, its influence on chaos is certainly an important key to 
understand the intrinsic mechanism of quantum chaos. In this paper, we only present the 
evidence of Pauh effect on chaos in a specific system. The most interesting and challenge 
questions are whether the dynamical Pauli effect can always enhance chaotic motions and 
what generic role does the Pauli principle play in nonlinear dynamics of quantum theory? 
We leave these questions for future investigations. 
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APPENDIX A: SP(6) COHERENT STATES 



To make this paper self-contained, we briefly list in this appendix the main results of 
the Sp(6) coherent states which has been discussed in Ref. |]T3[. The coherent states of the 
group USp(6) (the compact group of the Lie algebra Sp(6)) are defined by 



\v) = nvm 



(Al) 



where the generalized displacement operator T{ri), the coset representative of USp(6)/U(3), 
is given by 



fiv) = exp (^VooS++ E V2,D;^ - H-C. 



(A2) 



and |0) is the core vacuum (nearest closed shell), and -D^ the pair creation operators in 
the Sp(6) model. The matrix representation of T{ri) has the form: 



f (r/) Rep(f ) 



Y X 

-X+ Y 



(A3) 



where the matrices X and Y are the 3x3 subblock matrices of the coset representative of 
USp(6)/U(3). Their explicit forms are given by 



X 



22 



|X21 



1X21 



Xoo 



'Xoo + V^X 
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|X2-1 



(A4) 



X00 + V|X20 V 1X2-1 V3X2-2 

where xa^ is a function of {rjXfj,}. There are totally twelve parameters in this coherent 



state. Since the matrix T is unitary, the matrix X must be symmetric and the matrix Y is 
hermitian and is determined by X from the following relations: 



YY + XX^ 



YX = XY . 



(A5) 
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APPENDIX B: GEOMETRICAL MODEL VIA THE SP(6) MODEL 



In this appendix, the geometrical model and the Sp(6) model will be compared. Without 
any loss of generality, we define the nuclear surface variables as 

«2oW = /3{t) cos[7(t)] , 

«22(t) = ^/5(t)sin[7(t)], (Bl) 

where P{t) and 7(t) are the usual geometrical deformation parameters. Substituting Eq. ( PT]) 
into Eq. (0), it gives 



\Si\ = -aoo + V2P cos ^7 - -^^r^ , 



[B2) 



where i = z,x, y; = 0, kx = I, ky = 2 and 



The effective potential of the Sp(6) model is 

Veff = {H{7roo, 7r20, 7r22 = 0; ctoo, "20, "22)) , (B4) 

where tta^ is the conjugate momentum of axfi- When these variables are zero, the variables 
Si in Eq. (|T^ become real (for a proof, see Appendix C). Substituting the real variables Si 
of Eq. (p^ ) into Eqs. (^8| - p9l) , the effective potential can be expressed as: 



K//[G] = AP^ + BP^ + Caoop^ cos 37 + D + F . (B5) 

where the stiffness coefficients A, B, C, D and F are listed as follows: For G = SU(2), 

A = 7fi(fi/3 - l)/2 , 5 = -2r^(^]/3 - 1) , 

C = 2V2Vl{Vl/?,-1) , D = n{VL-n/2 + ?,n/{2VL))/Q, 

F=in/3r E S,S,[{l-S^){l-S])f^\ (B6) 

and for G = SU(3), 
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A = -7^](2^]- 1)/4 
c = -V2n{2n - 1) 



B = n{2n - 1) , 

D = 5n{2n-n)/{4n) 



F = 0. 



(B7) 



These results agree with that of Ref. [|T^]. The effective potential of the Sp(6) model also 
has the same form as the potential energy of the geometrical model in the intrinsic frame 



VefflP, 7] = IC,P^ + - J ^P' COS 37 + ■ 



(B8) 



In other words, the Sp(6) model in this picture contains all the properties of the nuclear 
geometrical model. 



APPENDIX C: DEFORMATION PHASE SPACE 

In this appendix, we will show that the surface variables, aoo? «2o and a'22 and their 
conjugate momentum ttoo, 1^20 and Tc'22 form the canonical variables, where 0:22 = "\/2ct22- To 
do so, we introduce the fourth type generating function 

It transforms the canonical variables {Vi, Qi, i = x,y, z} mentioned in Sec. Ill into another 
set of the variables {aA^(,7rA^}. The connection between these two sets of variables are 

3 P + 71"V2^"j' 

1 / , ^20 , /s ' \ 
3 P + 71^V2"^V' 

i (ttoo - v^vrao) , (C2) 

and = gives the same relation listed in Eq. (]TP|). The explicit connection between 
the variables Si and {ax^, t^\^) can be found by using Si = Vi exp(— i^). Hence {a\^, tta^} 
are the canonical variables. Eq. ( |C2|) and Si = Pj exp(— i-p^) imply that Si is real when 
T^Xfj. = 0. Thus we can easily obtain the effective potential discussed in Appendix B. 



Qy 



23 



REFERENCES 

[1] O. Bohigas, in "Chaos and Quantum Systems", Proc. NATO ASI Les Houches Summer 
School, Eds. M. J. Giannoni, A. Voros and J. Zinn- Justin, (Amsterdam, Elsevier, 1991). 

[2] M. Gutzwiller, "Chaos in Classical and Quantum Mechanics", ( Springer- Verlag, 1990). 

[3] W. M. Zhang and D. H. Feng, Phys. Rep. 252, 1-100 (1995). 

[4] W. M. Zhang, D. H. Feng, J. M. Yuan, Q. Pan, and T. Jjon, Phys. Rev. A 42, 3646 
(1990); W. M. Zhang and D. H. Feng, J. Mod. Phys. A8, 1417 (1993). 

[5] W. M. Zhang and M. T. Lee, "A simple picture of quantum suppression of chaos". 
Preprint of Institute of Physics, Academia Sinica, Taiwan. 

[6] A. K. Pattanayak and W. C. Schieve, Phys. Rev. Lett. 72, 2855 (1994). 

[7] B. Hu, B. Li and W. M. Zhang, "The quantum Frenkel- Kontorova Model, a squeezed 
state approach", submitted to Phys. Rev. Lett. (1997). 

[8] B. Hu, B. Li, J. Liu, and J. Zhou, "Quantum fluctuation: suppression or enhancement 
of chaos", submitted to Phys. Rev. Lett. (1997). 

[9] W. M. Zhang, D. H. Feng, and R. Gilmore, Rev. Mod. Phys. 62, 867 (1990). 

[10] W. M. Zhang, D. H. Feng and J. M. Yuan, Phys. Rev. A 42, 7125 (1990). 

[11] C.-L. Wu, D. H. Feng, X.-G. Chen, J.-Q. Chen, and M. W. Guidry, Phys. Rev. C 36, 
1157 (1987); C.-L. Wu, D. H. Feng, and M. W. Guidry Adv. Nulc. Phys. Vol. 21, 227 
(1995). 

[12] D, H. Feng, C.-L. Wu, M. W. Guidry, and Z. P. Li, Phys. Lett. B 205, 156 (1988). 

[13] A. Bohr and B. R. Mottelson, "Nuclear Structure", Vol. 2, (Benjamin, New York, 1975). 

[14] W. M. Zhang, C.-L. Wu, D. H. Feng, J. N. Ginocchio, and M. W. Guidry, Phys. Rev. 
C 38, 1475 (1988). 

24 



[15] W. M. Zhang, D. H. Feng, C.-L. Wu, and M. W. Guidry, J. Phys. G 15, L 115 (1989). 

[16] W. M. Zhang, C. C. Martens, D. H. Feng, and J. M. Yuan, Phys. Rev. Lett. 61, 2167 

(1988) ; W. M. Zhang, D. H. Feng, J. M. Yuan, and S. J. Wang, Phys. Rev. A 40, 438 

(1989) . 

[17] W. M. Zhang, and D. H. Feng, Phys. Rev. C 43, 1127 (1991). 

FIGURES 

FIG. 1. The Poincare surface of section [Q^ — 0) of the Sp(6) model with Q — 36, 
n = 10 and E = -3.0 MeV: (a) Gq = -0.1 MeV, 62 = 0. MeV, (b) Go = 0. MeV, 62 = -•! 
MeV, (c) Go = -0.01 MeV, 62 = --09 MeV. 

FIG. 2. The Poincare surface of section (Q^ — 0) of the Sp(6) model with energy 
increasing, where Q = 36, Gq = -.035 MeV, 62 = -.065 MeV: (a) n = 4, (b) n = 8, (c) 
n = 10, (d) n = 12, (e) n = 24, (f) n = 30. 

FIG. 3. A phase diagram of chaotic patterns in E — n plane in the Sp(6) model with 
Q = 36, where R means the regular region, W the weak-chaotic region, and S the strong- 
chaotic region. 
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